The origin of jerky dislocation motion in high-entropy alloys

Dislocations in single-phase concentrated random alloys, including high-entropy alloys (HEAs), repeatedly encounter pinning during glide, resulting in jerky dislocation motion. While solute-dislocation interaction is well understood in conventional alloys, the origin of individual pinning points in concentrated random alloys is a matter of debate. In this work, we investigate the origin of dislocation pinning in the CoCrFeMnNi HEA. In-situ transmission electron microscopy studies reveal wavy dislocation lines and a jagged glide motion under external loading, even though no segregation or clustering is found around Shockley partial dislocations. Atomistic simulations reproduce the jerky dislocation motion and link the repeated pinning to local fluctuations in the Peierls friction. We demonstrate that the density of high local Peierls friction is proportional to the critical stress required for dislocation glide and the dislocation mobility.

High-and medium-entropy alloys (HEAs and MEAs) are a new class of metallic materials that contain multiple elements at high concentrations and form solid solutions 1 in contrast to conventional alloys that typically consist of a single principal element with low concentrations of secondary elements. Since entropy is not always the decisive design parameter, they are a subclass of concentrated random alloys. Several HEAs outrival the mechanical properties of conventional alloys, some possess exclusive property combinations such as high strength and high ductility-even down to cryogenic temperatures [2][3][4][5][6] .
The origin of the high strength of HEAs is key for materials design but is still controversially discussed in the community [7][8][9][10][11] . There are, however, four different origins explored in the literature to date: i. Short-range ordering (SRO), i.e., a preferred ordering in the atomic-sized neighborhood, is well known to lead to a strength increase in conventional alloys 12 . This strengthening mechanism has been reported for the CoCrNi, a subsystem of the Cantor (CoCrFeMnNi) alloy 7,9,13 . Since trains of dislocations traveling on the same plane were observed in this system, it is argued that the first dislocation gliding on a plane destroys SRO within that plane and facilitates the glide of subsequent dislocations 14 . Yin and Curtin 8 , on the other hand, showed that the extraordinary strength of CoCrNi 2 can also be explained based on solid solution strengthening without resorting to SRO. ii. Indeed, dislocation pinning and the resulting strengthening are also observed in ideally random alloys 2 . Li et al. 7 reported that even in a perfectly random CoCrNi alloy some preferential atomic arrangements, in the case of Co and Cr, exist, and breaking these randomly appearing favorable bonds requires additional energy. iii. The unexpected wide range of stacking fault (SF) energies observed for HEAs could provide an alternative explanation for the pinning experienced by dislocations 15,16 . The resulting wavy dislocation lines have been observed for different HEAs in transmission electron microscopy (TEM) experiments [16][17][18] and atomistic simulations 7,19,20 . Li et al. 7 argued that the wide range of SFE causes peaks in the energy landscape where lattice friction becomes enhanced locally. Therefore, dislocation segments need to unzip from the local barriers while gliding. Such a mechanism would only exist in HEAs, but not in conventional dilute solid solutions as their local fluctuations in the SF energy remain comparatively small. iv. Lastly, the solid solution strengthening model developed by Varvenne et al. 21 predicts the flow stresses in an HEA solely based on misfit volumes and elastic properties. The dominating mechanism of solid solution strengthening in HEAs is intricate as it becomes impossible to define distinct matrix and solute atoms. A Labusch-type weak interaction 22 between solute atoms and dislocations is assumed in Varvenne's model. Here, the dislocation line interacts with the whole field of solutes within its elastic interaction range at the same time. This is opposite to the strongpinning model from Fleischer 23 where each solute corresponds to a site where the dislocation line is pinned.
In this study, we investigate the repeated local pinning along subsequent dislocations in the face-centered cubic (FCC) CoCrFeMnNi HEA by combining in-situ TEM deformation, atomic resolution analytical scanning transmission electron microscopy (STEM), and atomistic simulations. In-situ TEM deformation tests reveal wavy dislocation lines and a jagged dislocation motion under an applied tensile strain, although no indication of SRO is found by atomic resolution STEM energy-dispersive X-ray spectroscopy (EDS). Atomistic computer simulations reproduce the jerky dislocation motion with repeated pinning in random samples without SRO. The underlying mechanism is directly linked to local fluctuations in Peierls (force/friction) barriers, which lead to dislocation pinning points in concentrated random alloys. We further show their relation to the critical stress needed to move a dislocation line and the resulting mobility. While not a full solid solution strengthening model, the atomic fluctuations in frictional forces provide predictive design guidelines to tailor the dislocation pinning in HEAs.

Electron microscopy characterization of dislocation glide and pinning
We firstly analyze the structure and glide motion of dislocations in CoCrFeMnNi HEA by in-situ TEM tensile straining ( Fig. 1 and Supplementary Movie 1). Uniaxial tensile load is applied to the Cu grid which transfers stress to the TEM lamella prepared by focused ion beam (FIB) lift-out. The loading direction and the viewing direction are [113] and ½21 1, respectively (Supplementary Fig. 1) and three {111} slip planes are inclined with respect to the viewing direction. Due to the finite thickness of the TEM thin foil, which is around 100 nm, dislocation glide on the slip planes is constrained by the top and bottom surfaces. Figure 1a presents the representative TEM snapshots showing the glide motion of an array of Shockley partial dislocations which are separated by SFs appearing as dark bands. The dislocations bow out under the local resolved shear stress toward their glide direction. All slip system activated under tensile load. The TEM bright-field image was obtained at two beam diffraction conditions with strong excitation of g 022 beam. The SF bound by a pair of leading and trailing partial dislocations is visible with dark contrast. Difference images in the bottom row were prepared by using two TEM frames captured at denoted time intervals. Only a portion of the dislocation lines, marked by the red and yellow arrows, move within the time span. b Dislocation line of the leading partial dislocation (marked by "D" in a) extracted from consecutive TEM images. The dislocation lines show jerky glide motion; the dislocations exhibit a meandering shape as some segments get pinned during glide. c-e Series of atomistic simulation results showing the strong pinning of dislocations (red) during glide in the equimolar Cantor alloy at cryogenic (5 K) and ambient (300 K) temperatures and the weak pinning of dislocations in CoNi alloy at 300 K. images are processed to enhance the dislocation contrast (see section "Methods"). Tracing the evolution of the dislocation line shape during glide reveals that dislocations move through a kink-pair-like mechanism 24 . For example, only the upper part of the dislocation marked by a red arrow moves forward, then with some delay, the bottom part (yellow arrow) follows. The difference images in the right panel of Fig. 1a, highlighting the relative changes between two frames, show that only certain segments of dislocations move within a given time frame while others remain pinned.
This jerky motion can be seen more clearly by tracing the dislocation position over time (Fig. 1b). The positions are extracted using Gaussian and ridge detection image filters. As described above, a short segment of the dislocation line propagates first by forming a small hump, while the rest hardly moves. Most dislocations exhibit this type of localized pinning-unpinning along their lines constantly during glide, which could be the reason for the high shear stress required for dislocation glide in HEAs 17 . The effects of experimental artifacts, such as FIB damages or surface oxide, are expected to be marginal since these are present at the outer surface of the TEM sample with a small thickness of 5-10 nm 25,26 while the total thin foil thickness is around 100 nm. Therefore, the jerky motion in the middle of a dislocation, which is placed near the center of the TEM foil should not be affected by these surface artifacts.
In order to investigate whether the inherently meandering dislocation line shape seen in TEM can be reproduced via atomistic computer simulations, samples of comparable size are built. Each sample contains about 100 million atoms and features an inclined glide plane measuring 100 × 160 nm 2 (shown in Supplementary  Fig. 2a, c, d). To facilitate dislocation nucleation under tensile load a surface notch, aligned with this glide plane, is cut into the sample before loading the sample. Figure 1c-e show the first Shockley partial dislocations emitted and gliding away from this notch under tensile load. At low temperatures (5 K), strong pinning of the dislocation line in the HEA can be seen (Fig. 1c). These manifest as "V"-shaped dislocation line arrangements, where the line is pinned in one spot, while neighboring dislocation segments advance. These pinning points occur frequently, which results in a rough dislocation line. Thermal activation at ambient temperature renders dislocations apparently straighter (300 K Fig. 1d). The thermal energy lifts the dislocation over some of its weaker pinning sites, thus only the strong pinning points, which have a much lower density, persist active. For comparison, dislocations in the binary CoNi alloy were studied in an identical geometry (see Fig. 1e). There are almost no dislocation pinning points in this material, leading to straight dislocation segments during glide. The cusp in the middle of the dislocation stems from the sequence of initial nucleation events at the notch (cf. Supplementary Fig. 2d, e). Thus, the effect of strong dislocation pinning in HEAs and weaker dislocation pinning in binary alloys can be resolved. Moreover, the molecular dynamics (MD) simulations provide clear evidence for intrinsically wavy dislocations in the absence of SRO as observed, experimentally.
A comparison of experiment and simulation in Fig. 1 is only valid if both samples are in a similar random solid solution state. While it is clear from the sample preparation that there is no elemental clustering in the virtual sample, this question has not been answered definitively for the TEM sample. Lee et al. 17 previously performed atom probe tomography measurements on the samples shown in Fig. 1a and found no elemental ordering or clustering on the length scale of~1 nm. For further detailed analysis, we carried out atomic-resolution STEM EDS chemical mapping and quantitative analysis combined with EDS simulations  to determine the concentration of each atomic column around a Shockley partial dislocation. The high-angle annular dark field (HAADF) STEM image of the Shockley partial and its Bragg filtered analysis are given in Fig. 2a, b, confirming the a 0 /6[211] Burgers vector of the selected dislocation core, where a 0 is the lattice constant. With this quantitative evaluation of composition at each atomic column combined with EDS simulation, we could conclude the chemical homogeneity around the dislocation core without apparent clustering or segregation of specific element(s). There is fluctuation in per-column concentration though, which is, however, in line with the pristine bulk and negligibly small ( Supplementary Fig. 4). This result is also consistent with the recent high-resolution EDS experiments by Heczko et al. 27 , who report no elemental segregation to the strain field of Shockley partial dislocations in heat-treated Cantor alloy. A local descriptor for dislocation pinning in HEAs Figure 1 shows an inherent waviness of dislocation lines during glide in both experiment and simulation, indicating the existence of dislocation pinning sites. While we know that the simulated sample is in a perfectly random state without elemental clustering or short-range order, the experiments also show that there is no clustering or segregation around the dislocation core of Shockley partial dislocation.
Therefore, some statistically occurring atomic arrangements must be present, even in the random alloy, which pin the dislocation lines and locally hinder their glide. In the following, we investigate a single (dissociated) edge dislocation in the equimolar CoCrFeMnNi HEA to determine the nature of these nanoscale dislocation pinning points.
The established Peierls model of dislocation glide states that the stress required to move a dislocation is proportional to the maximum gradient along the energy path connecting two adjacent stable dislocation positions. For pure metals, this energy pathway (Peierls energy) may be approximated by a sinusoid and the so-called Peierls stress can be calculated analytically as a derivative of this energy. In dilute solid solutions, the local energy landscape is modified in the vicinity of each solute and therefore the local stresses required for dislocation glide changes 28,29 .
As the exact form of the Peierls energy landscape is usually unknown, the generalized stacking fault (GSF) energy surface 30,31 is often taken as a stand-in. The GSF energy surface maps the energy required to rigidly displace two densely packed crystal planes against each other. Supplementary Fig. 8 shows the GSF energy surface for a HEA sample. Here, the displacement and energy profiles corresponding to leading and trailing partial are indicated. As the GSF energy profile is obtained by averaging over the whole lattice plane, it is symmetric and smooth, showing no effect of the local chemical inhomogeneity. However, given the strong variance in local chemical arrangements, it is expected that the GSF energy varies on a local scale in HEAs.
It is widely accepted that the stable SF energy in HEAs, which marks a single point on the GSF curve strongly depends on the chemical environment 32,33 . Therefore, one would expect this dependence also for other points on the GSF landscape. From the locally fluctuating GSF landscape, one can directly conclude that the Peierls stress and thereby resistance against dislocation glide in HEAs becomes a spatially varying quantity. In the following, we will calculate this for each atom in the glide plane.
While the calculation of a local GSF energy landscape seems counterintuitive at first glance, it becomes a rather well-defined problem in the framework of classical atomistic simulations. LAMMPS 34 , by relying on interatomic potentials, is not only able to calculate the total energy for an arrangement of atoms but also assign a specific energy contribution to each atom. If we now perform a GSF energy calculation for the HEA system, we can store these per-atom energies and therefore obtain not only the global GSF energy landscape but also one resolved for each atom. While the general idea of this approach is similar to the stochastic Peierls-Nabarro model proposed in refs. 35,36, we do not attempt to find a closed form mathematical description of the local Peierls stress but instead measure the atomic Peierls force directly and confirm their effect on the dislocation pinning in a given alloy. Figure 3a shows the energy landscape for a single Co atom in the GSF plane. Compared to the globally averaged GSF energy landscape ( Supplementary Fig. 8) one can already see that the projected atomic energies are not symmetric anymore, as the local chemical environment changes after the initial displacement (red). Therefore, displacing this Co atom along the leading (red) and trailing (orange) partial dislocation displacement path does not result in the same energy barrier. Figure 3b shows the corresponding energy profiles along the two displacement paths. Dashed lines indicate the maximum gradient of the GSF curve, which corresponds to the local Peierls friction force. This value is proportional to the strength of this atom acting as a pinning point for the dislocation line and is denoted as F P . Note, this pinning point's strength is given in units of force (eV Å −1 ). It cannot be normalized to stress (eV Å −3 ) as the area per atom is not well defined in HEAs.
The next question is whether this approach actually works in describing the dislocation pinning on the atomic scale. To answer this, we insert a single edge dislocation into this pre-calculated GSF plane and apply a constant shear load to it. Figure 3c shows the leading partial dislocation of the dissociated edge dislocation gliding from left to right across the calculated pinning point landscape (the atomic pinning point strength F P is averaged on a 3 × 3 Å 2 grid and smeared out for better visibility). Movie S2 shows this sequence in motion. Here, the dislocation is pinned on the highest strength pinning point, i.e., the strongest, obstacle and bows out under the applied shear. Once the a, b GSF energy projected on a single atom in the glide plane. Red and orange arrows indicate the respective displacements connecting the local minima. The local Peierls barrier F P (in units of force) is defined as the maximum slope of this atomic energy curve. For reference, the GSF landscape averaged over the whole sample is given in Supplementary Fig. 8. c Representation of the local Peierls force in the glide plane (color corresponds to the local pinning force) with the dislocation (red) superimposed during glide. The per-atom F P values are averaged on a 3 Å-by-3 Å grid and smoothed using Gaussians for visual representation. It can be seen that the dislocation is pinned on the strongest obstacle, bends around it, and finally detaches. Note how the displacements caused by the passing dislocation change the pinning point landscape for subsequent dislocations. The mapped region spans 45 Å by 90 Å. d Superposition of different snapshots of the leading partial dislocation during glide. As snapshots are taken at constant time intervals, dislocation lines being closely together correspond to long dwell times, i.e., stronger localized pinning. The position of the strongest dislocation pinning points is marked by gray circles. Good agreement between the two can be seen. adjacent dislocation line segments have advanced a sufficient distance, the increase in line tension lifts the dislocation line across this strong obstacle and the dislocation continues its glide. After the dislocation has passed a given area the atoms are shuffled and the dislocation pinning point landscape changes for following dislocations.
Lastly, Fig. 3d shows the extracted position of the leading partial dislocation during glide over a 16 ps time interval. As the dislocation snapshots are taken every 250 fs, multiple closely spaced line segments correspond to the positions where the dislocation line is locally pinned while wider line spacing indicates faster dislocation glide. Superimposed are gray circles indicating positions of the highest local pinning point strength predicted from the atomic GSF calculations. Even by eye, a good match between proposed pinning points and regions of strong pinning is discernible. Obviously, there is no simple 1:1 correspondence between dislocation line shape and proposed pinning sites. This is the case as the "history" of the dislocation changes its current behavior. If the dislocation line was initially pinned and is therefore strongly bowed out, there is an additional force arising from the line tension to straighten out the dislocation line and shorten its overall length which can assist the crossing of pinning points along its way.
Given that our simulations have found an atomic scale descriptor for dislocation pinning points that is in line with the established theory of the Peierls model, we test this descriptor for a wide range of HEAs and search if an increased number of high F P lattice sites occurs.
Following the steps outlined in the previous section, we prepare a wide range of samples of different composition ranging from binary and ternary alloys to different HEAs. First, the atomic GSF landscape is calculated for all atoms in a given {111} plane. Figure 4a shows the calculated F P distribution g(F P ) for CoNi, CoCrNi, and CoCrFeMnNi. This comparison reveals that in CoNi almost all pinning points have a resistive force below 0.2 eV Å −1 . The addition of Cr leads to a substantial change in this spectrum. The interaction of Co and Cr causes the emergence of a new peak with F P values~0. 8  The three examples presented in Fig. 4a suggest the Co-Cr interaction to be responsible for strong dislocation pinning points. To confirm, we explore the mean concentration in the first and second nearest neighbor shell of the strongest (95th percentile F P,95 ) and the weakest (5th percentile F P,5 ) pinning points. The resulting concentrations, given in Fig. 4b, reveal that the concentration of Cr is increased around the strongest pinning points. Be reminded that the strong pinning points are exclusively Co atoms, confirming the interaction of these two species. The weak pinning points, on the other hand, mirror the bulk concentration (labeled "All") suggesting that there are distinct atomic environments responsible for the strongest dislocation pinning, even in a random alloy. However, the chemical concentrations of the atomic coordination acting as strong pinning sites show a great variance suggesting that there is not a single atomic arrangement responsible for the dislocation pinning.
Next, we investigate the effect of the pinning point spectrum on the critical force required to initiate dislocation glide as well as their mobility. One would expect dislocation pinning to be dominated by the fraction of strongest pinning points holding back the dislocation line. Therefore, we take F P,95 as a descriptor. For reference, this 95 percentile is marked in Fig. 4a. The dislocations in each sample are subjected to increasing shear until they start to glide. The critical applied force to initiate glide, F C , is determined as the force required to move the dislocation by 20 Å. These simulations are carried out at cryogenic temperatures to reduce thermal noise on the dislocation line shape. Note that this force is not necessarily identical for leading and trailing partial dislocations, especially for CoCrFeMnNi and CoCrNi which have negative SF energies, resulting in an additional driving force acting on the leading partial. Figure 4c shows F C over F P,95 for all samples. Blue symbols indicate results for the smaller samples with an initial dislocation line length l 0 of 20 unit cells, while red symbols indicate results from larger samples with an initial line length of 200 unit cells, i.e., lattice constants a 0 . A clear linear relationship between the pinning point density and the critical force can be seen within the small and large sample series confirming the expected correlation, where an increased density of strong pinning points leads to an increased force required to move a dislocation. Note that both sample sets are not on the same line, as the longer dislocation has a higher probability of having a locally weakly pinned segment where glide can start. This segment subsequently drags the other dislocation segments along. The mechanism is comparable to the weakest-link model proposed by Nöhring and Curtin 37 for cross-slip in FCC HEA, as suggested by Rizzardi et al. 38 .
We also calculate the dislocation mobility under a constant applied shear (Fig. 4d). Here, each sample is subjected to constant shear equivalent to the critical force required to move the leading partial dislocation by 10-30 Å (depending on the sample). In some samples, repeated pinning and depinning of the dislocation line were observed. Here, we isolated segments of constant velocity v and calculated the mobility M as, v = MF 28 , where F is the applied shear force. While the critical force was found to be proportional to F P,95 , the dislocation mobility should be inversely related to this value-the stronger the pinning points are, the slower and less mobile the dislocation will be. Therefore, we plot the extracted dislocation mobilities against F À1 P,95 . Again, a linear correlation between the dislocation mobility and the pinning point descriptor can be seen within each set of samples. Similar to the critical force, the weakest-link argument can be taken as an explanation for the differences between the two differently sized samples. Recent work by Sills et al. 39 showed a similar effect in stainless steels, where very long dislocation line lengths were required to converge the calculated dislocation mobility. However, comparing samples of identical initial dislocation line lengths should be valid.
To confirm that the correlations shown here are no coincidence, we also compare F C and M against other descriptors derived from the per-atom GSF curves. These can be found in Supplementary Figs. 12 and 13. Here we show the atomic stable and unstable stacking fault energies to F P, plotting their distribution following Fig. 4a and the resulting correlations akin to Fig. 4c, d. This confirms that F P is by far the best descriptor for F C and M.
In summary, we show that the dislocation pinning points identified in the Cantor alloy could also be found in different binary and ternary alloys. The critical force required for a dislocation to move correlates with the density of high strength pinning points identified from the atomic scale resolved GSF curves. Similarly, a greater number of strong pinning sites leads to overall reduced dislocation mobility in all investigated samples, independent of the material system or interatomic potential used. Investigation of the chemical environment around each pinning point reveals Co at its center and an increase in Cr in its surrounding, however, no single structural motif responsible for dislocation pinning has been confirmed. Given that the interaction of Co and Cr atoms needs to be harnessed, direct adjacency seems to be preferable.

Impact of the dislocation pinning points
We propose a descriptor for dislocation pinning sites in concentrated alloys. The descriptor is physically motivated and based on the atomic Peierls friction, which can be calculated as the maximum gradient of the per-atom GSF energy. We calculate the descriptor for a given pristine lattice plane, which is subsequently traversed by a dislocation allowing for direct spatial correlation of the established descriptor and observed dislocation pinning sites. Concentration analysis of the strongest dislocation pinning points reveals that they stem from the unique, nonlinear interaction of Co and Cr atoms.
The energy and force landscape felt by a dislocation is schematically shown in Fig. 5. In a pure metal, the energy change E encountered by a dislocation during glide is negligibly small. Therefore, the gradient of this energy and hence the force F required to move the dislocation is almost zero as well 40 . If a low concentration of solutes is added to the metal, forming a dilute solid solution, these solutes can either form local energy minima that trap the dislocation or high energy barriers blocking dislocation advancement. Both effects lead to steep energy gradients and therefore require high forces to be passed (Fig. 5a). In a concentrated HEA (Fig. 5b), on the other hand, the energy variance in the matrix is already higher hindering dislocation glide. However, as shown in Figs. 3 and 4, there are some environments that pin dislocations exceptionally strongly. These could stem from atomic arrangements where an energy minimum is located next to a barrier (as shown around d = 55 in Fig. 5b), leading to a steeper energy gradient. While these arrangements are theoretically also possible in dilute alloys, their probability is much lower given the lower density of secondary atoms. To confirm this schematic image, we also prepared linear profiles of the F P landscape in Co, Co 99 Cr 1 , and CoCrFeMnNi samples. These are shown in Supplementary Fig. 11 and confirm that there are no strong pinning points in the pure metal, some discrete strong pinning sites in the dilute CoCr solid solution and an abundance of pinning spots in the full HEA. The strong pinning spots in the CoCrFeMnNi alloy do not only show an increase in frequency but also higher peak F P values.
The identification of dislocation pinning points allows us to sample their density for different binary and ternary alloys and compare those to the CoCrFeMnNi HEA. As shown in Fig. 4a, neither Co nor Cr alone leads to a substantial shift in the F P distribution when added to Ni. Only the addition of both gives rise to a pronounced increase in pinning point density and strength. Further addition of Fe and Mn leads to no relevant changes in the F P spectrum and therefore only decreases the density of strong pinning points. The intuitive trends expected from the densities of strong pinning points, i.e., increasing resistance against dislocation glide and reduced mobility as the fraction of strong pinning points in a material increases are confirmed by our single dislocation shear simulations. Even though this is by no means a solid solution strengthening model, one might naively assume that an increased resistance against the onset of dislocation glide corresponds to an increase in the strength of a material. Under this assumption, our findings would be in line with experimental observations, where CoCrNi consistently shows higher hardness 41 and tensile yield strength 2,42 compared to other alloys from the same family, including the 5 component CoCrFeMnNi Cantor alloy. Additional simulations show that Ni is not even a necessary constituent, as a FCC random solid solution of CoCr has an even higher dislocation pinning than CoCrNi. This is in agreement with established literature, where FCC Co 63 Cr 29 alloys are commonly used in dental prosthetics due to their high yield strength and work hardening 43 .
Li et al. 7 reported the strengthening in regions of increased Co-Cr bonding by calculating the activation energy for dislocation motion in samples with different degrees of Co-Cr local ordering. Their calculation method requires the use of dedicated smaller samples, which means that they cannot calculate the energy barrier for long-range dislocation glide. Our approach, on the other hand, allows the determination of the atomic pinning point strength in the glide plane of the dislocation. Therefore, a direct spatial correlation between the atomic descriptor and pinned dislocation line can be made. Antillon et al. 44 also proposed a method of calculating the local energy barrier against dislocation glide. Their method, however, relies on energy measurements as the dislocation passes through the material in a MD simulation. In contrast, the proposed descriptor can be calculated from the pristine material, i.e., without an embedded dislocation, under static conditions. Therefore, this descriptor has the possibility of being used for future materials design, as it does not require costly MD simulations of many hundred-thousands of atoms to identify potential chemical compositions of strong dislocation pinning.

Transferability of our results
The simulation results outlined in this manuscript were all obtained at cryogenic temperatures. This suppresses the thermally assisted passing of dislocation pinning points. Within the framework of transitionstate theory, the energy barrier (not the pinning forces considered in Figs. 3 and 4) associated with a pinning point becomes important to determine whether a dislocation is pinned on a given pinning site. Therefore, the pinning points found here will only pin a dislocation, if the applied force is below the pinning point strength and the associated energy barrier is greater than the thermal energy. Assuming that the energy landscape is unchanged by both temperature and the presence of the dislocation, we can take the per-atom unstable stacking fault energy as a surrogate for the local energy barrier. Considering our data, we can see that the unstable stacking fault energy shows broad distributions in these alloys: CoNi 73.2 ± 16.7 meV, CoCrNi 100.4 ± 32.9 meV, and CoCrFeMnNi 74.7 ± 47.7 meV, with overall relatively low energy barriers suggesting that thermal activation will become important at intermediate temperatures.
We propose that the dissociated edge dislocation studied here is representative of FCC materials. In general, the calculated F P values are intrinsically insensitive to the dislocation character, as they are calculated based on the pristine lattice. This also extends to the actual pinning point-dislocation interaction. Our descriptor only accounts for the rigid displacement of the atoms during the passage of the dislocation. This means that even though the shear direction (in relation to the line direction) is different for screw and edge dislocations, the shear process is comparable. Hence, there is no obvious dependence of our descriptor on the dislocation character. However, the displacement paths that need to be calculated to obtain F P depends on the shear direction, i.e., the Burgers vector. This is an obvious difference from established solid solution strengthening models, e.g., ref. 21, where the elastic interaction between solute and dislocation depends on their stress fields and therefore on the dislocation character. Based on these arguments, we also expect the coupling of Co and Cr, shown in Fig. 4b, to be independent of the dislocation character.
Dislocations in body-centered cubic (BCC) metals glide differently than dislocations in FCC metals. Usually, the Peierls barrier in BCC metals is much higher requiring a kink-pair mechanism to advance the dislocation 28 . Based on this mechanism, the displacement of the atoms is most likely different from the rigid shear displacements covered by our F P descriptor. Therefore, it is most likely not directly applicable. However, if one finds the rigid shear "path" or trajectories of the atoms involved in such a process in the BCC matrix, the resulting force barriers could be calculated in a similar way. The glide of screw dislocations through a BCC HEA matrix has recently been studied in detail by Wang et al. 45 , who use nudged elastic band calculations to understand the energy landscape of the gliding dislocation through the BCC matrix. An increased glide resistance can be seen around the energy walls and valleys. b shows E and F for concentrated HEAs. Here, the variance in E and F is higher in the matrix. However, a combination of adjacent local minima and maxima (d = 55) leads to a much higher variation in local pinning forces and higher maxima.
Our study reveals the atomic origin of dislocation pinning points in multi-principal element alloys thereby answering a long-standing question in the materials science and engineering community. We find that no distinct SRO is required for strong dislocation pinning in HEAs as local pinning of the dislocation line is facilitated by a locally increased slope of the GSF curve. This scenario agrees with the established Peierls model but reveals that HEA and specific alloys are susceptible to exceptional Peierls peak heights. In the alloy samples investigated, these changes are highly non-linear and the interaction of Co and Cr leads to the highest density of strong dislocation pinning points. This understanding of the local Peierls stress landscape and the resulting pinning forces acting on dislocations gives access to a different parameter that can be tuned in the design of alloys with increased resistance against the onset of dislocation glide and overall reduced dislocation mobility.

Fabrication of CoCrFeMnNi single crystal
A single crystal CoCrFeMnNi HEA was produced by using the Czochralski method from an equiatomic melt of high-purity elements. The composition of Fe, Co, Cr, Mn, and Ni is 24,24,20,12,20 at.%, respectively, measured by both STEM-EDS and atom probe tomography 17,46 . In situ TEM straining samples were prepared by focused ion beam (FIB) site-specific lift-out method 47 that utilizes a 30 and 5 kV Ga ion beam (JIB-4610F, JEOL). A thin lamella with a size of several micrometers was glued on a custom-made Cu grid by Pt deposition.

In situ TEM tensile experiments
In situ tensile straining tests were carried out in an aberrationcorrected TEM operated at 300 kV (JEM-ARM300CF, JEOL) equipped with a straining holder (Model 654, Gatan). In situ TEM movies were recorded with a high-speed camera at a frame rate of 50 frames per second and 2k resolution (Oneview, Gatan). The displacement of the whole~9 mm large sample support was controlled by a step motor whose resolution was 1 μm. The straining was interrupted intermittently to stabilize the stage and observe dislocation glide motion. No quantitative load-displacement data can be obtained during the deformation.
The contrast of dislocations in TEM images and videos is enhanced by image processing. Firstly, the background noise from the FIB damage is removed by using direct subtraction of the last frame after aligning all frames using the cross-correlation method. Then, a Gaussian filter is used to blur the remaining noise. To extract dislocation positions as shown in Fig. 1b, a ridge detection filter is used after applying an additional Gaussian filter.

STEM EDS experiments
The atomic-resolution STEM EDS chemical mapping was carried out on a JEM-ARM 200F (JEOL) equipped with a spherical aberration corrector (ASCOR, CEOS) and energy dispersive X-ray spectrometer (JED-EDS, JEOL). Multi-frame EDS maps were acquired for a half-hour by using dual-type EDS detectors (the effective X-ray detection area of a 100 mm 2 for each detector) with a large effective solid angle (~1.2 sr) and a highly focused electron probe (~1.2 Å) at the electron dose rate of 6.3 × 10 9 e nm −2 s −1 . The resulting elemental maps were obtained by the multiple frame summation up to 2300 frames with 256 × 256 pixels resolution and the acquisition time of 10 μs per pixel (~20 min in maximum as a total acquisition time). Multi-frame EDS maps were averaged by correcting the frame-to-frame displacement caused by sample drift during acquisition. The averaged X-ray count of each element was converted to the composition by calibrating the k-factor, which accounts for the different X-ray yields of each element. The background noise floor in each composition map was suppressed by averaging every neighboring 5 pixels.

Atomistic simulations
The equimolar CoCrFeMnNi HEA and its subsystems were all simulated based on the 2NN MEAM 48 interatomic potential by Choi et al. 49 , the CoCrCuFeNi 20,50 and the Co 30 Fe 17 Ni 37 Ti 17 19,51 HEA are based on the Farkas and Caro 50 or Zhou et al. 51 EAM potentials 52,53 , respectively, to exclude erroneous correlations from a given parametrization.
All simulations were run in LAMMPS 54 . Samples were prepared using ATOMSK 55 , and post-processing was based on algorithms implemented in OVITO 56 , and accelerated by PARALLEL 57 .

Virtual TEM lamella
The large-scale sample shown in Fig. 1 measurements had initial dimensions of 101 × 152 × 75 nm 3 of pristine FCC lattice oriented 1 23 Â Ã , 1 2 1 Â Ã , [412] along x, y, and z, respectively. This gave approximately 10 8 lattice sites, filled with either CoCrFeMnNi or CoNi alloy. A surface notch with an inclination of 28.13°was cut into the y-z plane. This surface notch was aligned with a {111} glide plane in the sample. To mimic the experimental conditions, the x-and z-directions were taken as open surface boundaries, while the y-direction was set to periodic boundary conditions. This sample geometry is schematically shown in Supplementary Fig. 2a.
After equilibration for 50ps at 5K or 300K in the NVT ensemble, the sample was subjected to uniaxial strain along y up to a strain of 0.06 using an engineering strain rate of 10 8 s −1 58, and a constant time step length of 1 fs. Dislocation lines were extracted every 5 ps during straining using the dislocation extraction algorithm DXA 59,60 .

Shear of an edge dislocation
To cut a dislocation into the samples, a half-plane was inserted into an FCC crystal resulting in a misfit edge dislocation following the geometry shown in Supplementary Fig. 2b 61 . The resulting dislocation had a Burgers vector b along 1 10 Â Ã (corresponding to the x-direction) with a line direction of 11 2 Â Ã (corresponding to the y-direction). Two different sample dimensions were considered for the isolated dislocation simulations. The smaller sample contained 477600 atoms with periodic boundary conditions along 1 10 Â Ã (x) and 11 2 Â Ã (y) directions and an open boundary along the [111] (z) direction. This sample was built from 200 × 20 × 20 unit cells in each of those directions, respectively. The larger sample had the same crystallographic orientations and boundary conditions, but the total size was increased to 400 × 200 × 24 unit cells leading to~11.5 × 10 6 atoms. In both cases, the lattice sites were filled with the desired atomic species to create a random solid solution.
All samples were annealed for 50 ps at 2 K using a 1 fs time step. During this annealing and the subsequent shear simulations the surface layers were thermostatted using a Langevin thermostat 62 , while the atoms in the center of the sample were integrated in the NVE ensemble.
After equilibration, a ramping shear force was applied to each atom in the outer surface layer. This force pointed in the x-direction (c.f. Supplementary Fig. 2b) and was ramped from 0 to 0.208 eV Å −1 over 200 ps. From these shear ramp simulations, the critical force to initiate dislocation glide FC was determined.
Once this critical force was established, a second constant force simulation was started. Here, the applied shear force was first ramped from 0 to FC with a rate of 10 −6 eV Å −1 ps −1 . Afterward, it was held for 100 ps.
The dislocation position was extracted using DXA 59,60 implemented in OVITO 56 .

GSF surface calculation
The global and atomic GSF curves were calculated for the two differently sized samples presented in the previous section before dislocation insertion (Supplementary Fig. 2b).
For a GSF calculation, the sample is split into two crystallites. One is above the <111> GSF plane, while the other one is located below this plane. Applying the conventional simulation methodology prevents the atomic relaxation during rigid displacement of the two crystallites within the GSF plane 31 , however, the relaxation of the atoms normal to the GSF remains unconstrained. One problem that arose, applying this methodology, was that the atoms in the disordered alloy relaxed into different intermediate states. While this effect averages out in the global GSF calculation, the resulting atomic energy landscapes are not smooth but include discontinuities when an atom relaxes into such a metastable state. To avoid these effects, we further constrained the out-of-plane atomic relaxation during the calculations. During displacement, the forces acting on the atoms along the 1 10 Â Ã (x) and 11 2 Â Ã (y) directions 31,63 were set to 0 while the forces normal to the GSF plane were averaged for all atoms in each crystallite. This means that the pristine FCC lattice was preserved within each crystallite. The crystallites, however, were still able to change their separation distance to reduce the system's energy. Supplementary Fig. 9 shows the difference between this rigid relaxation method compared to the conventional atomic relaxation scheme (relaxed).
Even though this constrained relaxation approach does not resolve the effects of the intrinsic lattice distortions found in HEAs directly, it does show their effect implicitly. A comparison of Fig. 3a and Supplementary Fig. 8 reveals that the local energy minima in the GSF landscapes do not have to be in the same positions for the global GSF landscape and the ones resolved for each atom, even though both are extracted from the same simulation. The spatial variations in the equilibrium positions for each atom are caused by local fluctuations in chemical environments. In the subsequent calculation of the energy gradients these minima are connected, which implicitly accounts for the lattice distortions.
To determine the GSF landscape, a 3.6 × 2.1 Å 2 area was sampled with a grid spacing of 0.1 Å. At each grid point, the potential energy was minimized using the latest implementation of the FIRE algorithm 64 to an energy tolerance of 10 −12 eV. The resulting two-dimensional energy landscape for each atom was interpolated using cubic splines and the three local minima were found using the BFGS algorithm (both implemented in SCIPY 65 . The maximum energy gradient was calculated along a linear path connecting these three minima. Figure 3 shows an example of the resulting atomic energy landscape and calculated energy gradient.

Data availability
All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Information. Additional data related to this paper may be requested from the authors.